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ABSTRACT 

Planets are formed in disks around young stars. With an age of ~ 10 Myr, TW Hya is one of the nearest T Tauri stars that is still surrounded by a 
relatively massive disk. In addition a large number of molecules has been found in the TW Hya disk, making TW Hya the perfect test case in a large 
survey of disks with Herschel-PACS to directly study their gaseous component. We aim to constrain the gas and dust mass of the circumstellar 
disk around TW Hya. We observed the fine-structure lines of [O I] and [C II] as part of the Open-time large program GASPS. We complement this 
with continuum data and ground-based '^ CO 3-2 and ''CO 3-2 observations. We simultaneously model the continuum and the line fluxes with 
the 3D Monte-Carlo code MCFOST and the thermo-chemical code ProDiMo to derive the gas and dust masses. We detect the [O I] line at 63 yum. 
The other lines that were observed, [01] at 145 yum and [CII] at 157 yum, are not detected. No extended emission has been found. Preliminary 
modeling of the photometric and line data assuming ['-CO]/[''CO]=69 suggests a dust mass for grains with radius < 1 mm of ~ 1.9 x lO"'' Mq 
(total solid mass of 3 x 10"' Mq) and a gas mass of (0.5-5) x 10"' Mg. The gas-to-dust mass may be lower than the standard interstellar value of 
100. 

Key words. Circumstellar disks 
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1. Introduction 

Planets are formed in the disks that surround a large fraction of 
T Tauri stars. Knowledge of the gas mass available at different 
disk ages is essential to constrain giant planet formation mod- 
els. Most studies estimate the dust mass from millimeter contin- 
uum emission and assume the gas mass is a factor of 100 times 
larger. This conversion factor has been calibrated for the inter- 
stellar medium but is likely not valid for disks, especially those 
that are evolving toward debris disks or where most of the gas 
has accreted onto the planetary atmosphere. Disk gas mass esti- 
mates derived from observations of '"CO and optically thinner 
'^CO emission are at least a factor of 10 lower than the mass de- 
rived from dust observations assuming the interstellar medium 
conversion factor The discrepancy has been ascribed to CO pho- 
todissociation at disk atmospher e and freeze-out onto cold dus t 
grains in the disk midplane (e.g.. lOi et al.ll2004tlThiet al.ll2001h . 
An alternative explanation is that the CO abundance is not dif- 
ferent and the gas in disks has been deplete d. 



servations of all the major gas-phase carbon and oxygen-bearing 
species, we can more precisely constrain the disk gas mass. 



The PACS instrument jPoglitsch & al.l 2010h on-board the 
Herschel Space Telescope (iPilbratt & al.l2010l) makes it possible 
to observe lines from species that result from the photodissocia- 
tion of CO (atomic oxygen and singly ionized carbon). With ob- 



* Herschel is an ES A space observatory with science instruments pro- 
vided by Principal Investigator consortia. It is open for proposals for 
observing time from the worldwide astronomical community. 



At a distance of ~ 56 pc ( IWichmann et al.lll998h . TW Hya is 
one of t he nearest classical T Tauri st ars with an estimated age of 
10 Myr dBarrado Y Navascuesll2006l) . Its proximity allows us to 
attain an order of magnitude higher mass sensitivity than objects 
in the Taurus molecular cloud. Fits to the spectral energy distri- 
bution (SED) provide an estimate of the gas disk mass of 6 x 
10"^ M q after applying a conversion factor of ~ 75 (ICalvet et al.l 
l2002h . This large disk mass at this advanc ed age is surprisin g 
as the median disk lifetime is only 2-3 Myr dHaisch et al.ll200lT) . 
TW Hya is considered a transition object wi th an op tically thin 
inner cavity and an optically thick outer disk JCalvet et al. 2002t 
iRatzka et al.ll2007l) . The fit to the SED also suggests that grains 
have grown to at least ~ 1 cm. 



The star TW Hya was observed as a Science Demonstration 
Project object and is part of the Herschel-GASPS program 
(iDent & GASPS teami l2010l) . Herschel observations of the 
disk around the H erbig Ae star HD 169 142 are presented by 
'Meeus et al.l (l2010l) . In this letter we use fine-structure lines in 
addition to continuum data and CO (sub)millimeter lines to di- 
rectly constrain the gas mass and compare it to the dust mass 
derived from fits to the SED. 
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Table 2. Disk parameters for the modeling. 



63.0 63.1 63.2 

Wavelength [|J,m] 



63.3 



63.4 



Fig. 1. Herschel-PACS spectrum centred around the OI 63 yum 
hne on the upper-left panel. 

Table 1. Lines observed by Herschel-PACS. The errors and up- 
per limits are 3 cr. The calibration error adds an extra ~ 40% 
uncertainty. The CO data also have uncertainties of 30%. 



Line 


Cont. flux 


Obs. GH08 


M08 




(Jy) 


(10-'** Wm-2) 




01 3Pi ^'Pj 


2.99 ±0.14 


36.5 ± 12.1 124-161 


412 


01 3Po ^3pj 


7.00 ± 0.05 


< 5.5 25-41 


11 


CIl2p3/,^2p,/2 


8.79 + 0.08 


< 6.6 0.8-12 


0.06 


CO 3-2 


n.a. 


0.43 0.3-0.6 


n.a. 


I'CO 3-2 


n.a. 


4.4 X 10-' n.a. 


n.a. 







Fixed parameters 










Stellar mass 
Stellar luminosity 
Effective temperature 
Solid material mass density 
Inner radius 
Outer radius 
ISM UV field 


M.(Mo) 
L.(Lo) 

7-eff(K) 

Pdu,si(g cm ') 
R,u(AV) 
RoUAV) 
X 


Inner cavity 

0.6 

0.23 

4000 

3.5 

0.25 

4 

1.0 


Outer ring 

4 
196 


a viscosity parameter 
Turbulent velocity 
Disk inclination 




a 

v,urb(kms-') 

i 






0.0 
0.05 

7 




CO isotopologue ratio 




['^COM'^CO] 






69 






MCFOST best fit parameters 








Column density index 
Reference scale height 
Reference radius 
Flaring index 




6 

Wo(AU) 
7 




2.0 
100 
0.6 


1 


10.0 
100 
1.12 


Minimum grain size 
Maximum grain size 
Dust size distribution index 
Dust mass (a <1 mm) 
Solid mass 


«min(/Jm) 

aniax(cm) 

P 

A^du,sl(MG) 
Wsol,d(Mo) 


1.2 

2.C 


X 10 
X 10 


X 10 
10 

3.4 

-9 
-8 


_2 

1.9 X lO-'' 
3.0 X 10-3 




ProDiMo parameter range 








Disk gas mass 
UV excess 

Fraction of PAHs w.r.t. 
Cosmic ray flux 


ISM 


Fu\ 
/pah 
f(s-') 




3 X 10-*-0.3 

0.018 

0.01,0.1 

(1.7-17) X lO-" 



2. Observations and results 

We obtained photometry in the "blue" (70 /im) and "red" 
(160 lum) band of the PACS camera by doing mini scan maps 
with a scan speed of 20" and a scan length of 2' (obsid 
1342187342). The total duration of this map was 731 sec, with 
an on-source time of 146 seconds. The results are 3.90 + 0.02 Jy 
and 7.38 + 0.04 in the blue and red band respectively and have 
an absolute accuracy estimated to be 5% for the blue channel and 
10% for the the red channel. These values agree very well with 
the observed IRAS flux densities and also with the continuum 
flux densities measured with the PACS spectrometer (Table [1]). 
We also used the PACS spectrometer to target the [OI] line at 
63 fiva in line scan mode, and the [OI] and [CII] lines at 145 
and 158yL(m, respectively in range scan mode (obsid 1342187127 
PacsLineSpec and obsid 1342187238 PacsRangeSpec). Only the 
[OI] line at 63 //m was detected and we report upper limits for 
the other two lines; see Table[T] The absolute accuracy of PACS 
spectroscopy is currently estimated to be about 40%, but is ex- 
pected to improve in the future. Figure [T] shows the spectrum 
centered at the position of the 01 line at 63 fim of the central 
pixel. 



3. Modeling and discussion 

As there is no evidence for an outflow from TW Hya, we as- 
sume that all the fluxes a rise from the circ umstellar disk (see 
also the discussion in iMathews et alJ 120101) . The interpretation 
of the observation s with the DENT grid of models is detailed in 
iPinte et aP (l2010l) . We performed a more detailed analysis here. 



We first augmented the Herschel photometric data with con- 
tinuum measurements from the literature. We also retrieved 
and reduced archival SCUBA data for TW Hya obtained dur- 
ing two nights with very good sub-millimeter transmission 
(Fv(450 //m)= 4.25 + 0.85 Jy and 7^^(850 //m)= 1.38 + 0.14 
Jy). The disk around TW Hya has an internal cavity from up 
to 4 AU where the gas and dust density are very low. Most 
of the mass is located in the external ring. The inner (/?{„) and 
outer radius (.Rout) of the external ring are well constrained by 
imaging studies and are fixed at 4 AU and 200 AU resp ectively 
dRoberge et al.ll2005HOi et ^ "2004'. "H ughes et al.ll2007l) . We fit- 
ted the SE P with the 3D M onte-Carlo radiative transfer code 
MCFOST ( iPinte et al.l |2006|) . We chose to restrict to a para- 
metric disk model for this letter. The disk has a radial density 
profile with index e. The flaring is characterized by an open- 
ing angle Hq at a given radius /?ief and a flaring index y so 
that the gas scale-height is given by // = HoiR/RrefY. The 
low continuum flux in the 30-100 jum region suggests that the 
outer disk flaring is weak. Amorphous olivine grains were used 
(iDorschner et al.l [1995) with a power-law size-distribution de- 
fined by a minimum radius «„,!„, maximum radius flmax, and 
power-law index p. The dust size-distribution and mass are well 
constrained by the continuum emission at long wavelengths. The 
fit to the long-wavelength photometric points including the new 
Herschel-PACS data is shown in Fig. |2] and the disk parame- 
ters constrained by the fit are listed in Table |2] The inferred 
dust mass in grains with radius < 1 mm is Mdust=l-9 x 10"^ 
Mq and the total mass in solids (pebbles) up to Omax = 10 cm 
is Msoiid=3 X lO-"* Mq. However, the fit fails to account for 
the flux at ~25 /im, which may stem from our assumption of 
a unique temperature for grains of all sizes. The flux around 
20-30 fiva is strongly inclination-dependent because we adopted 
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a sharp density change between the inner cavity and the outer 
ring at 4 AU. Solids as large as 10 cm in radius a re needed to 
account for the observed 7 mm and 3.6 cm flux dWilner et alj 
I2OOOI) . The small grains in the TW Hya disk account for 6% of 
the total solid mass. We also estimated a mass in small grains 
of 3 X 10""* Mq assuming that the emission in the millimeter is 
optically thin, an average dust temperature of 20 K, and grain 
opacity Ky - 2.0(v/vo)~^ g cm~^ wh ere vq - 230.769 GHz and 
P - 0.6 dBeckwifli & Sargenliri991h . The two estimates of dust 
mass (with radius a < 1mm) are consistent within a factor 2 with 
each other. The visibility amplitudes generated by the m odels are 
consi stent with the observed amplitudes at 345 GHz by dOi et alj 
12004 . 

For the line o bservations we augmented the Herschel data 
wifli SMA CO 3-2 dOi et al.l2004 and JCMT '^^CO 3-2 observa- 
tions dThi et al.l2004l) . Following the characterization of the disk 
structure from the SED, we ran three series of models with the 
th ermo-chemical cod e ProDiMo (a detailed description is given 
in lWoitke et alJ2009h . In ProDiMo species abundances are com- 
puted at steady-state from the gas, and dust temperature as well 
as the local UV field for the photodissociation reactions. A con- 
stant isotopologue ratio ['"'CO]/['~CO] of 69 is assumed. The gas 
kinetic temperature is computed by balancing heating and cool- 
ing processes. Line profiles are computed by non-LTE radiative 
transfer within ProDiMo. The disk is assumed to be passively 
heated. The disk turbulent velocity and incl ination are wefl con- 
strained by millimeter interferometric data dOi et al.l|2004l) . The 
outer disk is irradiated by direct and scattered stellar photons 
as well as by interstellar UV photons. The free parameters of 
the gas simulations are the disk gas mass Mgas (between 3 x 
10""^ and 0.3 M©), the fraction of polycyclic aromatic hydrocar- 
bons (PAHs) in the disk with respect to the interstellar abun- 
dance /pah, and the cosmic ray flux C, (=1.7 x 10"'^ s"' in the 
ISM). Observations show that PAHs are depleted by at least a 
factor of 10 (fpAH - 0.1) in d isks with respect to the interstellar 
abundance dGeers et al.ll2006l) . Because the gas is mostly heated 
by photoelectrons ejected from PAH, the PAH abundance is the 
main free parameter that controls the gas temperature. The three 
series of models correspond to three possible states: disks with 
a very low PAH abundance (/pah - 0.01), disks with a typical 
PAH abundance (/pah = 0.1), and X-ray irradiated disks with 
a low PAH abundance (/pah = 0.01) but ten times the standard 
cosmic ray flux {C-\.l x 10~'^ s~') to mim ic the influence of 
strong X-ray emission dBruderer et al.ll2009l) . The model results 
are plotted in Fig. [3] The density, dust and gas temperature struc- 
ture are shown for a typical disk in the appendix. The results 
from series 3 are within 10% of the values of series 2, suggest- 
ing that X-ray does not influence the line fluxes that are emitted 
at radii beyond a few AU. In panels a and b we can see that the 
01 63 jum and 145 //m flux increases with the disk gas mass. 
The 01 63 //m line is optically thick while the OI 145 jum line 
is optically thin for all models. Both lines arise mostly in a ring 
between 4 and 10 AU and thus probe the gas mass up to 10- 
20AU with 10-20% contribution from the inner cavity (panel f). 
In panel c the CII flux first starts to increase with higher disk 
gas mass but then plummets for disk gas masses greater than 
10"^ Mq. The CII line is optically thin and the flux increases 
with radius. As the disk becomes more massive, more carbon is 
converted into CO and the disk becomes cooler. Ionized carbon 
is excited in gas at ~ 100 K. The CO 3-2 flux increases with 
increasing disk gas mass although the emission line is highly 
optically thick with r >100 (panel d). CO 3-2 emission comes 
from the outer disk {R > 50 AU). Finally, panel e illustrates the 
use of the line emission ratio between two isotopologues ('-CO 
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Fig. 2. Fit to the SED generated by ProDiMo using the param- 
eters from MCF OST. The input Phoenix st ellar spectrum plot- 
ted i n red is fromlBrott & Hauschildtl d2005h . lUE (UV) data are 
fromlValentietalJ dlQOj. The 2MASS J,H, K, IRAS, and Spitzer- 
MIPS photo metry are arch i val da ta. The Spitzer-IRS spectrum is 
published byHatekaeTaD (l2007h . The Herschel-PACS data are 
plotted in filled green tria ngles. The average UBVRI p hotomet- 
ric points are published bv lRucinski & Krautterl(ll983h . The 800 
jjm. and 1.1mm data points (inverted filled blue triangles) are 
taken fromlWeintraub et al. (198^The 3.4mm point (open blue 
square) is from lWilner et al.l (12003) while the 7mm an d 3.6 cm 
points (filled blue square) are from lWilner et al.l (l2000h . 



and '-^CO) to constrain column densities or masses. The flux dif- 
ference between the two isotopologues shrinks with increasing 
disk gas mass. The observed '^CO/'^'CO 3-2 ratio is consistent 
with a very low-mass disk. The CII and CO lines probe the outer 
disk mass (panel f). All together, the observations constrain the 
disk gas mass between 5 x 10"^ and 5 x 10""^ M©. The modeling 
suggests that the disk of TW Hya has a gas-to-dust mass ratio of 
2.6-26, around a factor of 10 lower than the interstellar value. If 
we compare the gas mass to the total mass in solids (ie including 
solids with radii up to 10 c m), the gas-to-solid ratio is 0.17-1.7. 
iGorti & HoUenbachI d2008i GH08) included X-ray and UV heat- 
ing in modeling the disk of TW Hya with a gas mass of 0.03 Mq 
but noticed that X-ray weakly influences the fluxes. Their model 
overestimates the two OI line fluxes (Table [T] and Fig.[3]i. Disk 
models with X-ray hea ting only also predict too strong OI fluxes 
dMeiierink et al.ll2008l M08) for their model with Lx=2 x 10^° 
erg s ' scaled to d-56 pc. At ~ 10 Myr, TW Hya is one of the 
oldest classical T Tauri stars. The outer dust disk is very long- 
lived, while the inner disk contains little amount of material. The 
gas may have a shorter lifetime than the dust due to photoevap- 
oration or the small grains result from collisions between the 
large grains. TW Hya is one of the strongest X-ray active T Tauri 
stars ( Raassen 2009), which ma y result in a high gas photoevap- 
oration rate dOwen et al]|20IOl) as evidenced in t he blu eshifted 
[Ne II] emission observed bv lPascucci & SterzikI d2009h . 
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Fig. 3. Three series of model results compared to observations. The blue boxes enclose the model outputs for disk gas mass between 
5 X 10""* Mq and 5 x 10"'' Mq. Panel a shows the predictions and observation of the Ol 63 jiva line. The 3cr uncertainty range is 
plotted as dashed lines. Panel b and c show the predicted fluxes and the 3o- upper limits for the Ol 145 yum and Cll lines. The two 
lower panels (d and e) are the comparison between observations and model outputs for '^CO 3-2 emission and the '^CO/'^CO 3-2 
ratio. Panel f shows the normalized cumulative fluxes for a 10"^ M© model (series 1). The diamonds (o /?out=174 AU model, o 
^out=120 AU model) show the predictions for TW Hya from GH08. 



4. Conclusion 

The Herschel-PACS spectral observations were used to constrain 
the gas disk mass surrounding the 10 Myr T Tauri star TW Hya. 
We estimate the gas mass to be (0.5-5) x 10""' M© compared 
to the dust mass (amax < 1mm) of 1.9 x 10""* M©. The gas-to- 
dust mass ratio is ~2.6-26, lower than the standard interstellar 
value of 100. The ratio gas-to-total-mass in solids is ~0. 17-1.7. 
Although the disk is still massive, a significant fraction of the 
primordial gas has akeady disappeared. A large fraction of the 
primordial gas may have been evaporated due to the strong X- 
ray flux from TW Hya. TW Hya is the first example where the 
disk gas mass around a transitional T Tauri star can be deter- 
mined accurately and directly from gas phase lines. However, 
more detailed modeling that includes X-ray physics and '^CO 
photochemistry is needed to confirm the low gas mass. 
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Appendix A: Density and temperature structure 



We show in Fig. lA.lllA.2l and IA.3l the density, dust temeprature 
, and gas temperature profile respectively for a disk model with 
Mgas=2.9x lO-^M0,/pAH=O.l,and^= 
parameters are given in Table |2] 
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Fig. A.l. Density profile , 
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Fig. A.2. Dust temperature profile. 
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Fig. A.3. Gas temperature profile. The contour of Av=l is shown 
in white. 



